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Recent Developments in Heavy- Quarkonium Phenomenology 

Hee Sok Chung**, Jungil Lee*** and Chaehyun Yu** 
'^Department of Physics, Korea University, Seoul 136-701, Korea 

We review recent developments in heavy-quarkonium phenomenology within the nonrelativistic QCD factor- 
ization approach. Main issues we consider in this work include the polarization of prompt J/ip at the Fermilab 
Tevatron and the large relativistic and QCD corrections to double-charmonium production at the B factories. 
We also consider inclusive charm production in bottomonia decays. 



1. Introduction 

For more than a decade, the nonrelativistic 
quantum chromodynamics (NRQCD) factoriza- 
tion approach [1] has been employed to predict 
various measurables in heavy-quarkonium phe- 
nomenology. In this approach, rates for produc- 
tion or decay of a heavy quarkonium that in- 
volve creation or annihilation of a heavy-quark- 
antiquark {QQ) pair are expressed in linear com- 
binations of NRQCD matrix elements, which 
reflect the long-distance nature of the heavy 
quarkonium. Corresponding short-distance co- 
efficients are insensitive to the long-distance in- 
teraction and calculable perturbatively. As the 
first application of the factorization theorem for 
the decays, the light-hadronic and electromag- 
netic decay rates of the P-wave heavy quarkonia 
were factorized [2] . Explicit infrared singularities 
appearing in order-a^ color-singlet-modcl predic- 
tions for the light-hadronic decays of P-wave 
heavy quarkonia [3,4] became to be identified as 
those of the next-to- leading-order (NLO) correc- 
tions to the long-distance color-octet NRQCD 
matrix elements that mix with the color-singlet 
NRQCD matrix elements at this order. 

Soon after the resolution of the infrared 
problem in the P-wave quarkonium decay, 
'0(25') anomaly was reported: the transverse- 
momentum {pt) distributions for prompt '4>{'2jS) 
and J/il) measured by the CDF Collaboration at 
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the Fermilab Tevatron were greater than theoret- 
ical predictions of the color-singlet model by or- 
ders of magnitude at large pr [5,6]. Based on the 
idea of gluon- fragmentation dominance of quarko- 
nium production at large pr [7], the puzzle was 
resolved by introducing the color-octet 5'-wave 
spin-triplet [QQb,(^ Sx)] contribution to the gluon 
fragmentation into an S'-wave spin-triplet heavy 
quarkonium [8] . The dominance of the color-octet 
mechanism in large-pT prompt ipi'^S) production 
led to the prediction that the prompt ^{2S) must 
be transversely polarized at large pr [9,10]. Un- 
fortunately, the predictions of the polarizations 
of prompt 'il}{2S) [11,12] and prompt J/ip [13,14] 
were disfavored by the CDF measurement at Run 
I of the Fermilab Tevatron [15]. For more details, 
we refer the readers to Refs. [16,17]. This problem 
has been a serious challenge to heavy-quarkonium 
theory. 

Another interesting problem in quarkonium 
phenomenology that has not been resolved un- 
til very recently is regarding the cross section 
for the exclusive J/ip -\- rjc production in e^e~ 
annihilation at B factories. The cross section 
measured by the Belle Collaboration [18] was 
greater than theoretical predictions at leading or- 
der (LO) in Us and v by an order of magni- 
tude [19,20]. The BABAR Collaboration also 
observed a large cross section [21] and subse- 
quent studies [22,23] of the Belle Collaboration 
revealed that the large discrepancy between the- 
ory and data is not originated from the contam- 
ination of other signals [24-26]. On the other 
hand, there were arguments that an appropri- 
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ate choice of light-cone distribution amplitude 
may resolve the problem [27,28] , which were con- 
fronted with an opposing argument that, after 
separating the short-distance contribution from 
the long-distance light-cone distribution ampli- 
tude, the light-cone method yields a result equiv- 
alent to that of NRQCD [29]. Very recently, 
large relativistic corrections [30,31] in combina- 
tion with one-loop QCD corrections [32,33] filled 
the gap within errors [34]. 

In this work we review recent developments 
in heavy-quarkonium phenomenology within the 
nonrelativistic QCD factorization approach, fo- 
cusing on the following three selected topics. Is- 
sues regarding the polarization of prompt J/tp at 
the Fermilab Tevatron and future opportunities 
to resolve the problem at the CERN LHC are pre- 
sented in Sec. 2. The discrepancy between the LO 
predictions and the empirical cross sections for 
e+e~ J/ip+r]c measured at the B factories has 
been a long-standing puzzle. In Sec. 3 we briefly 
review the problem and describe a recently de- 
veloped method of resumming a class of relativis- 
tic corrections to all orders in the heavy-quark 
velocity v, which played a crucial role in resolv- 
ing the e+e~ — > J/ip + rjc problem in conjunction 
with the rclative-ordcr-a,, corrections. As the last 
topic, inclusive charm productions in bottomonia 
decay being analyzed by the CLEO Collaboration 
and corresponding theoretical predictions are re- 
viewed in Sec. 4 which is followed by conclusions 
in Sec. 5. 

2. Polcirization of prompt J /if) 

The discrepancy between the theoretical pre- 
diction and the CDF measurement at Run I of 
the Fermilab Tevatron for the prompt S-wave 
spin-triplet charmonium produced at large pT 
was embarrassing. Gluon fragmentation has the 
kinematic enhancement factor of p^/m^^^ which 
dominates over the suppression factor a'^{pT) for 
the coupling compared to the gluon-fusion dia- 
grams at sufficiently large pr [7]. In compari- 
son with the gluon fragmentation into the color- 
singlet QQ pair [QQi{^Si)], the QQsi^Si) contri- 
bution has a less severe suppression factor of order 

which accounts for the spin-non-flipping long- 



distance transition from the color-octet to the 
color-singlet state, while the suppression factor 
a^(pT) in the color-singlet fragmentation chan- 
nel is for the short-distance coefficient [8] . If the 
gluon fragmentation into the QQsi'^Si) pair is 
the dominant contribution to the prompt J/ip 
and ip{2S) production at large px, then the in- 
variant mass for the fragmenting virtual gluon 
is negligible compared to the typical scale pr 
of the short-distance process. Hence, longitudi- 
nally polarized virtual-gluon contribution must 
be suppressed compared to the transverse con- 
tribution with a multiplicative scaling factor of 
m^j/^/Pr [9,10]. Because the long-distance tran- 
sition from QQs.iC'Si) to QQii^Si) is dominated 
by spin-non-flipping chromoelectric dipole transi- 
tions, prompt '0(25') must be transversely polar- 
ized at large pt- This is also supported by a lat- 
tice computation of the color-octet NRQCD ma- 
trix element involving the spin-flipping interac- 
tion for the decay [35]. On the contrary, the pre- 
dictions for the ijj{2S) polarization [11 13] were 
disfavored by the CDF data at both Run I and 
Run II of the Fermilab Tevatron [15,36]. 

In a similar manner as iJj{2S), one can expect 
that the prompt J/ip, which has much larger 
statistics than tp{2S), produced at large pr must 
be transversely polarized. However, prompt J/4' 
suffers from feeddowns from higher resonances 
such as ip{2S) and Xcj- At the Tevatron, the 
feeddowns from tp{2S) and XcJ are about 10% 
and 30%, respectively. The contribution from 
XcJ dilutes the polarization and that from iJj{2S) 
is dominantly spin preserving [13,14]. Although 
feeddown from higher resonances makes the pre- 
diction for the prompt J/0 less dramatic than 
'ij){2S), the tendency of rising transverse polar- 
ization as Pt increases is the same as that for 
ijj{2S) [13,14]. Again, this expectation was con- 
fronted with the CDF data at both Run I and 
Run II of the Fermilab Tevatron [15,36]. One 
may suspect any flaws in the theoretical predic- 
tions. The prediction for the prompt J/ijj polar- 
ization in Ref. [13,14] was at LO in and v. 
One may expect that the inclusion of real gluon 
emission and higher-order contributions in v that 
flip spin may dilute the polarization even if they 



Recent Developments in Heavy Quarkonium Phenomenology 



3 



are suppressed based the velocity-scaling rules of 
NRQCD. One may also suspect that the domi- 
nance of gluon fragmentation has not yet been 
achieved in the pt range accessible at the Fermi- 
lab Tevatron. At the same time, it is fair to notice 
that the discrepancy between the Run I and Run 
II data of the CDF Collaboration is not clearly 
understood, yet [15,36]. 

Theoretical prediction for the T{nS) produced 
at the Fermilab Tevatron was not dramatic as the 
S-wSkve spin-triplet charmonium state because 
the available range of px does not reach the re- 
gion where gluon fragmentation dominates [37]. 
Theoretical uncertainties were large because of 
the large errors in various NRQCD matrix ele- 
ments of bottomonia [38]. The CDF measure- 
ment of the T{nS) polarization [39] was con- 
sistent with the prediction within errors [37]. 
Around the time when the prompt J/tl) polar- 
ization measurement of the CDF Collaboration 
was announced, the E866 Collaboration reported 
that T(IS') is almost unpolarized while the sam- 
ples for T{2S) + T{3S) are almost transversely 
polarized [40] in their fixc;(i-targc;t experiment. 
This shows that the polarization of T(IS'), that 
has many channels of feeddown from higher res- 
onances through P-wavc states, is diluted com- 
pared to T(2S') and T(3S') states. Similar ten- 
dency has been observed by the DO Collaboration 
at Rim II of the Tevatron [41]. The DO result for 
T(IS') shows that T(IS') is longitudinally polar- 
ized at most of Pt regions below about 15 GeV. 
This may be due to the feeddown effect through 
P-wave resonances. However, transverse contri- 
bution increases continuously from the maximally 
longitudinal region around 2-5 GcV as pT in- 
creases. The DO result for T(2S'), which has less 
feeddown effect than T(15) shows that the rise 
of transverse polarization is faster and stronger 
than that of T(IS'). This indicates that polar- 
ization analysis for T(n5) at the CERN LHC, 
where accessible px may be large enough to ap- 
proach gluon-fragmentation-dominance region for 
T{nS) production, may reveal transversely polar- 
ized spin-triplet S-wavc bottomonia. In compar- 
ison with the charmonium counterparts, T{nS) 
is more nonrelativistic. As a result, it will suffer 
less from the contamination of longitudinal con- 



tributions of higher orders in as and v. 

We list recent theoretical attempts that have 
been made to improve the NRQCD predictions 
for the inclusive heavy-quarkonium production at 
hadron colliders. The QCD NLO corrections to 
the color-singlet channel to the gluon-fusion dia- 
grams, gg — » QQi{^Si) + g, that contributes to 
the ^Si quarkonium production at hadron collid- 
ers were recently computed [42]. It was reported 
that the NLO corrections are large and mod- 
ify the expectations for the total cross sections 
and the distribution in the transverse momen- 
tum. This is a signal that the NRQCD matrix el- 
ement for the long-distance process QQsi^Si) 
J/ijj + X may be smaller than that estimated be- 
fore, resulting in less dramatic increase of trans- 
verse polarization of prompt J/tp. In Ref. [43], 
the contribution of the following specific chan- 
nels, gg — > QQi{^Si) + QQ, where Q is c or b, 
were calculated. The authors reported that, in 
those channels, fragmentation approximation un- 
derestimates the cross section in the kinematical 
region accessible at the Fermilab Tevatron. In ad- 
dition, the polarization of both J/'iJ; and T pro- 
duced through above channels are unpolarized. In 
Ref. [44], the authors update the theoretical pre- 
dictions for hadroproduction of direct T. They 
reported that the QCD NLO corrections to the 
color-singlet contribution significantly enhances 
the theoretical prediction for the cross section at 
large pt and substantially affect the polarization 
of the T [44]. Because the measured polariza- 
tions are for the inclusive T that contain various 
feeddowns, one cannot make a direct comparison 
with the available data for the inclusive T{nS). 
Further theoretical predictions for partial contri- 
butions to the inclusive charmonium production 
at hadron colliders are also available [45-47]. 

3. Relativistic corrections to quarkonium 

processes 

In general, because of clean initial states, pro- 
cesses for quarkonium production in e+e~ anni- 
hilation may have less theoretical uncertainties 
than hadroproduction. Especially, the exclusive 
double- charmonium production from e+e~ anni- 
hilation into a virtual photon should be described 
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relatively with ease because the color-singlet con- 
tribution in NRQCD is absolutely the dominant 
production mechanism. For the exclusive process 
e+e~ — »■ J/ip + r]c, only a single long-distance 
color-singlet NRQCD matrix element, {Oi)j/^, 
involves, once we impose the approximate heavy- 
quark spin symmetry of NRQCD. {Oi)j/^ is the 
most accurately known NRQCD matrix element, 
which is determined from T\J/0 c+e^]. 

However, the cross section first measured by 
the Belle Collaboration, which is about 33 fb [18], 
was greater than the NRQCD predictions [19,20] 
at LO in both and v by an order of mag- 
nitude. Later the Belle Collaboration updated 
the value as 25.6 ± 2.8 ± 3.4 fb [22] and the 
BABAR Collaboration also reported a large cross 
section of 17.6 ± 2.8 ± 2.1 fb [21], where the val- 
ues include the branching fraction of the sys- 
tem recoiling against the J/tp into two or more 
charged tracks. The theoretical predictions were 
3.78±1.26 fl) [19] and 5.5 fli [20], where the differ- 
ence depends on the choice of input parameters. 
The prediction in Ref. [19] also includes QED con- 
tributions that make a constructive interference 
of about 30% of the total cross section for the 
QCD process. As was described earlier, proposals 
that the signal may include final states other than 
J/ip + r]c [24-26] were disfavored by an updated 
analysis by the Belle Collaboration [22,23]. And 
the light-cone method reproduced the NRQCD 
prediction [29]. This puzzle has been a serious 
challenge to quarkonium phenomenology. 

When the first theoretical prediction was made, 
it was reported that there may be huge relativis- 
tic corrections to the LO prediction [19]. Unfor- 
tunately, there was no reliable way available to es- 
timate the NRQCD matrix element involving rel- 
ativistic corrections, which is power- ultraviolet- 
divergent and needs subtraction [19]. This was 
an obstacle to evaluating relativistic corrections 
to the process. As the next step, QCD NLO cor- 
rections to the e~^e~ J/iJj + rjc process were 
calculated. The K factor for the QCD NLO pre- 
diction to that of LO in both ag and v was about 
a factor of 2 [32,33]. In spite of a large enhance- 
ment from the QCD NLO corrections, large dis- 
crepancy between theory and experiments still re- 
mained. 



The only sector that had not been explored 
thoroughly was relativistic corrections to the pro- 
cess. The first prediction of the K factor for 
the relativistic corrections was 2.0l}°f^ [19]. The 
main source for the huge uncertainty was from 
the relative-order-u^ NRQCD matrix element. In 
principle, once the quarkonium color-singlet wave 
function is known by using a non-perturbative 
method such as lattice NRQCD [48,49], one can 
compute the order-t;^ NRQCD matrix element, 
where power-ultraviolet-divergent contributions 
to the matrix element must be subtracted with 
dimensional regularization to make it consistent 
with the standard NRQCD factorization formal- 
ism. Because of this large subtraction, even the 
sign of the matrix element was uncertain by the 
time the prediction was made [19]. 

As a breakthrough, a new method to determine 
the color-singlet S'-wave NRQCD matrix elements 
was introduced [30]. Employing the Cornell po- 
tential [50], which is a linear combination of the 
Coulomb potential and the linear potential, the 
authors computed the quarkonium color-singlet 
wave function [30,51]. Relative strength of the 
Coulomb potential compared to the linear poten- 
tial was determined by making use of the lattice 
string tension [30]. Furthermore, they derived the 
relation, (q^")// = (g^)^ [30], which is the gen- 
eralization of the Gremm-Kapustin relation [52], 
where (q^") is the ratio of the relative-order- 
v"^^ color-singlet matrix element to {0\)h [30]. 
The generalized Gremm-Kapustin relation made 
it possible for one to resum a class of relativis- 
tic corrections of color-singlet contributions to 
all orders in v [30]. The resummation method 
was applied to improve the NRQCD factoriza- 
tion formula for the electromagnetic decay of the 
J/ip. They determined the parameters of the 
Cornell potential model and the normalization of 
the quarkonium wave function by making use of 
the radial mass splittings of the 5-wave states 
and electromagnetic decay widths of the S'-wave 
states. As a result, the LO and order-t;^ NRQCD 
matrix elements of the 5-wave quarkonium H 
were determined [30]. In an updated analysis, 
they reported {Oi)j/^ = O.UOto fZ GeV^ and 
{q'^)j/^ = 0.441 GeV2 [53]. This result cor- 
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responds to {v'^)j/ij, = 0.225^0.088' which is con- 
sistent with the naive estimate « 0.3 within 
errors. This is the most rohablc determination 
of the color-singlet NRQCD matrix elements for 
J/ip available so far [53]. 

The method of resummation [30] has been ap- 
plied to resum a class of the relativistic correc- 
tions to the e+e" Jf-ip + rjc [31] in combination 
with the QCD NLO corrections in Ref. [32]. In 
the following updated study in Ref. [34] , the au- 
thors reported the prediction a[e'^e~ J/ip + 
ric] = 17.6lg;3 which agrees with the experi- 
ments within errors. For more detailed descrip- 
tions of estimating theoretical uncertainties, we 
refer the readers to Ref. [34]. Therefore, in spite 
of large error bars, it seems fair to conclude that 
the long-standing puzzle between theory and ex- 
periments for the e+e~ — > J ftp + r}c has been 
resolved [34]. 

Independently of the theoretical analysis in 
Ref. [34], a fixed-order relativistic corrections of 
relative-order to the cross section for e+e" — *■ 
J/ip + rjc is given in Ref. [55], where the non- 
perturbative factors {Oi)j/^, and {q^) j/^, are de- 
termined by making use of the measured values 
for V[J/tl) e+e-], T[ric 77], and T[J/tp 
?>g\ with the assumption {Oi) j/,ij{q^) = 
{Oi)n^{q^)nc- They obtained the cross sec- 
tion 20.04 fb with = 0.573 GeV^ and 
{q^)j/^, = 0.202 GeV^ [55]. The values for the 
nonperturbative factors yields {v'^)j/-^ ~ 0.09. It 
is smaller by about a factor of 2 than that ob- 
tained in Ref. [53], which is consistent with the 
naive estimate v"^ ~ 0.3. We notice that the short- 
distance coefficient of r[J/V' 3g] at relative- 
order is very large. This restricts the fit in 
Ref. [55] to have a small value for {v'^)j/^ [53]. 
Further applications to other exclusive double- 
charmonium production processes are also avail- 
able in Refs. [56,57]. 

4. Inclusive charm production in bottomo- 
nia decays 

As described in Sec. 1, the earliest calculations 
of the widths of P-wave quarkonium states us- 
ing perturbative QCD were plagued with infrared 
divergences [3,4]. In 1992, Bodwin, Braaten, 



and Lepage showed that the infrared divergences 
could be absorbed into the probability for the 
QQ pair to be at the same point in a color- 
octet state [2]. Based on the NRQCD factoriza- 
tion formalism, one can now carry out rigorous 
calculations of inclusive charm production from 
Xhj decays. On the experimental side, Xbj{^P) 
and XbjC^P) have been discovered. The only 
properties of these states that have been mea- 
sured thus far are their masses and their radia- 
tive branching fractions into T{nS). The total 
widths of the Xbj{nP) states have not been mea- 
sured. Recent runs of the CLEO experiment at 
the T(25) and T(3S') resonances have provided 
new data on the X6j(l-P) and Xbj{2P) states. 
The B-factory experiments BABAR and Belle 
can study the Xbj {nP) states by using data sam- 
ples of T(2S') and T(3S') provided by initial-state 
radiation. The Belle experiment has also accu- 
mulated data by running directly on the T(3S') 
state. Therefore, it is necessary for theorists to 
present quantitative predictions for the P-wave 
bottomonium states that can be compared with 
forthcoming data. 

The widths of all four states in a P-wave mul- 
tiplet can be calculated by using the NRQCD 
factorization formula, once the two nonperturba- 
tive factors (Ci)x6 and {Os)xb have been deter- 
mined. In Ref. [58], the authors computed the 
inclusive charm production rate in the Xbj decay 
XbJ —^c + X. They also provided the momentum 
distributions of charmed hadrons produced in the 
process by making use of a parametrization for 
the charm fragmentation function fit to the data 
from e+e~ collisions [59]. The predictions are de- 
pendent on the ratio ps = m'^{Os) Xb / i^i) Xb that 
can be determined phenomenologically. Recently, 
the CLEO Collaboration measured the total pro- 
duction rates for the charm in Xb.j decays and 
the branching fractions for Xbj(lP, 2P) D^X 
with puo > 2.5 GeV [60]. The CLEO measure- 
ment indeed determined p^{\P) = 0.160^qo47 
and P8(2P) = 0.074lJ^;°J° [60]. The joint fit for 
IP and 2P gives ps = O.OSetJ^ I^?^ [60]. 

As shown in Ref. [58], the branching fraction 
Br[x6J c + X] must be particularly large for 
the J = 1 state. The reason is that the Xbi ccg 
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is one of the LO contributions Xbi ~* QQ9 to 
hadronic decay of Xbi > while the mode Xbj — * ccg 
is suppressed by order as compared to the lead- 
ing contribution XbJ — * 99 for hadronic decays 
of J = and 2 states. Note that Xbi 99 
is forbidden by Yang's theorem. The measure- 
ment of the CLEO Collaboration also confirmed 
the large branching fractions of about 25% for 
both X6i(l-P) and Xbi(2-P) [60], which is consis- 
tent with the prediction in Ref. [58] within errors. 

Theoretical studies on the inclusive charm pro- 
duction has been extended to the 5-wave bot- 
tomonium states in Refs. [61-64]. 

5. Conclusions 

The discrepancy between the theoretical pre- 
diction and the CDF measurements at Run I 
and Run II of the Fermilab Tevatron for the 
prompt S-wave spin-triplet charmonium pro- 
duced at large pr is embarrassing. To make the 
problem more complicated, there is a big discrep- 
ancy between the CDF data of Run I and those 
of Run II that has not been clearly understood, 
yet. In order to find any other contributions that 
may modify the theoretical prediction toward ex- 
perimental results, efforts have been made in var- 
ious manners. A few sources which contribute 
to the longitudinal polarization for the prompt 
J/tp and ip{2S) are found. However, the com- 
plete analysis which includes all of these contri- 
butions and includes all the feeddowns that are 
contained in the data is not available, yet. In 
order to make a reliable prediction for the po- 
larization with these new contributions, more ac- 
curate determination of various NRQCD matrix 
elements for the quarkonium production must be 
completed, which may not be an easy task. 

The puzzle of the cross section for e+e^ —> 
J/tp + Vc at the B factories has spurred rapid 
developments in heavy- quarkonium phenomenol- 
ogy. It led the first reliable estimate of the 
relative-order-t!^ color-singlet NRQCD matrix el- 
ement for the J/ip, that was uncertain even in 
the sign. A new method to rcsum relativistic; cor- 
rections to a class of color-singlet contributions 
has been developed to determine the color-singlet 
NRQCD matrix elements for the 5-wave heavy 



quarkonia more accurately. A recent computation 
of order-as corrections to the quarkonium electro- 
magnetic current to all orders in v revealed that 
the velocity expansion for approximate J/ ip oper- 
ator matrix elements converges very rapidly [54] . 
This supports the reliability of the resummation 
method. Being applied to the e+e~ — » J/tp + Vc 
process, the resummation method for the rela- 
tivistic corrections has resolved the puzzle in com- 
bination with large QCD corrections. The puzzle 
also motivated the proofs of the factorization the- 
orems for exclusive two-body charmonium pro- 
duction in B-meson decay and e+e~ annihilation 
to all orders in perturbation theory in quantum 
chromo dynamics [65]. 

Resolution of the infrared-singularity problem 
in P-wave heavy-quarkonium decay into hadrons 
was the first successful theoretical application of 
the NRQCD factorization. However, the process 
has never been thoroughly investigated experi- 
mentally. A new measurement by the CLEO Col- 
laboration of the inclusive charm production rate 
in P-wave bottomonium decay has provided use- 
ful information. The measurement helped us de- 
termine the ratio for the color-octet NRQCD ma- 
trix element to the color-singlet matrix element 
for the P-wave states accurately. On-going anal- 
ysis by the CLEO Collaboration on the charm- 
hadron momentum distribution may reveal the 
color-octet mechanism in P-wave bottomonium 
decay more clearly. 
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